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Abstract. We present the results of sub-mm, mm (850/xm, 450 pm and 1250/xm) and radio (1.4 and 4.8 GHz) 
continuum observations of a sample of 27 K -selected Extremely Red Objects, or EROs, (14 of which form a 
complete sample with K < 20 and I — K > 5) aimed at detecting dusty starbursts, deriving the fraction of 
UltraLuminous Infrared Galaxies (ULIGs) in ERO samples, and constraining their redshifts using the radio-FIR 
correlation. One ERO was tentatively detected at 1250/rm and two were detected at 1.4 GHz, one of which has 
a less secure identification as an ERO counterpart. Limits on their redshifts and their star forming properties 
are derived and discussed. We stacked the observations of the undetected objects at 850/rm, 1250/rm and 4.8 
GHz in order to search for possible statistical emission from the ERO population as a whole, but no significant 
detections were derived either for the whole sample or as a function of the average NIR colours. These results 
strongly suggest that the dominant population of EROs with K < 20 is not comprised of ULIGs like HR 10, but 
is probably made of radio-quiet ellipticals and weaker starburst galaxies with L< 10 12 Lq and SFR< 100 Mg 
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- Radio continuum: galaxies - Submillimeter 


1. Introduction 

The existence of a population of extragalactic objects with 
extremely red infrared-optical colours has been known for 
a number of years (for a recent review, see Cimatti 2000). 
These objects were initially discovered mainly in near- 
infrared surveys of blind fields and quasar fields and are 
very faint or invisible in the optical bands (Elston et al. 
1988; McCarthy et al. 1992; Hu & Ridgway 1994). These 
extremely red objects (EROs) are defined as those which 
have R — K colours > 5 and these tend to have K mag¬ 
nitudes fainter than ~18. Since the time of their discov¬ 
ery, the nature of this population has remained a puzzle. 
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Based on their NIR and optical photometric data, two 
broad classes of models are consistent with the observed 
red colours : (a) high redshift starbursts, red because of 
severe dust extinction. From the required extinction and 
simple dust models, these galaxies could be normal star- 
burst galaxies or even high-z counterparts of local Ultra 
Luminous Infrared Galaxies (ULIGs) and (b) old passively 
evolving ellipticals at redshifts greater than about one. 
The red colours of the EROs would then be explainable 
by a large R-correction and an absence of ongoing star 
formation. If the EROs belong to the former class, then 
they would be dominant sites of star formation and would 
be important in determining the star formation history of 
the universe (Cimatti et al. 1998a). On the other hand, 
if they belong to the latter class, then the volume den- 
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sity of these objects as a function of redshift would pose 
strong constraints on the models for the formation of el¬ 
liptical galaxies, which range from monolithic collapse to 
dark matter dominated hierarchical structure formation 
scenarios (Daddi et al. 2000b and references therein). 

Over the last few years, it has been possible to study 
in detail a handful of EROs which are bright enough to 
yield reliable spectra and have multi-wavelength contin¬ 
uum data. These studies were able to determine the na¬ 
ture of these EROs and there are now examples known for 
both starbursts (Cimatti et al. 1998a; Cimatti et al. 1999; 
Smail et al. 1999a; Smail et al. 1999b; Gear et al. 2001) 
as well as for old ellipticals (Spinrad et al. 1997; Cimatti 
et al. 1999; Soifer et al. 1999; Liu et al. 2000) among the 
ERO population. Studies indicate though, that not more 
than ~30 % of EROs are starbursts (see Sect. 5 for de¬ 
tails). Recently, independent wide-held surveys for EROs 
have been conducted (Daddi et al. 2000a; McCarthy et 
al. 2000; Thompson et al. 1999) which have shown that 
these objects are strongly clustered in the sky (see also 
Chapman et al. 2000 and Yan et al. 2000). The surface 
density of these EROs after correcting for clustering, as¬ 
suming that these are passively evolving ellipticals, has 
been shown to be consistent with pure luminosity evolu¬ 
tion with a formation redshift greater than 2.5 (Daddi et 
al. 2000b). 

HR 10 is one of the reddest EROs and is quite bright 
A'=18.42, Graham and Dey 1996). 
detected 850 /mi and 1250 /tm 
emission from this object (see also Dey et al. 1999). They 
derived a star formation rate of several hundred M 0 yr -1 
and an FIR luminosity in excess of of 10 12 L 0 , thus show¬ 
ing that HR 10 is an ULIG at its redshift of 1.44 (Graham 
and Dey 1996). At the time of the observations presented 
in this paper, HR 10 was the only ERO with detected 
submm emission and also the only ERO with a known 
redshift. Hence it was thought possible that a majority 
of EROs would be similar to HR 10 and would have ob¬ 
servable submm and mm emission (Cimatti et al. 1998b, 
Andreani et al. 1999). Therefore we began a search for mm 
and submm continuum emission from other EROs with an 
aim to detect extreme starburst galaxies. Further, in order 
to constrain the redshifts of these objects and also to un¬ 
derstand their nature, we also decided to search for radio 
continuum emission from these EROs. 

The radio and the FIR luminosities of nearby galax¬ 
ies (z < 0.4) which are dominated by star formation are 
known to be highly correlated (see Condon 1992 and ref¬ 
erences therein), whereas E and SO galaxies are known 
to be more radio bright than a star-forming galaxy of 
similar FIR luminosity (Walsh et al. 1989). The radio 
(mainly synchrotron and some amount of free-free) emis¬ 
sion and the FIR (due to dust) emission have different 
spectral indices. Hence, assuming the local radio-FIR cor¬ 
relation holds at high redshift, the observed ratio of radio 
to FIR emission strengths can be used to determine the 
redshift of a star-forming galaxy (the redshift determina¬ 


tion method and the associated error estimation are de¬ 
veloped in Carilli and Yun (1999, 2000) and Blain 1999). 

So, if on the one hand, the ERO population consists 
primarily of starbursts, this method could be used to de¬ 
termine the nature of EROs and also to constrain the red- 
shifts of these objects. If, on the other hand, the EROs 
are ellipticals, their IR colours would be used to constrain 
their properties. Since most EROs are too faint to obtain 
redshifts even with lOm-class telescopes or even to obtain 
accurate photometry over the entire optical-IR range, such 
complementary diagnostics become important in under¬ 
standing these objects. 

We first describe the ERO sample, the observations 
and their results in Sects. 2 and 3. In Sect. 4, we derive 
statistical properties of the sample and estimate the frac¬ 
tion of starbursts and ellipticals with radio emission in 
Sect. 5. Those EROs with radio or mm detections are dis¬ 
cussed in detail and their properties are derived in Sect. 
6. The cosmology adopted in this paper is a flat universe 
with Ho=70 km s -1 Mpc -1 and all results are calculated 
for both Ha=0.7 and Da=0.0. The spectral index is as¬ 
sumed to be —0.7 throughout this paper for calculating 
the A'-correction for the radio emission. 

2. Sample and Observational details 

A sample of 27 EROs was selected from the literature 
and also from other samples being studied by the authors. 
These objects were identified from a variety of observa¬ 
tions : deep NIR surveys of random fields, galaxy cluster 
fields and quasar fields (see references in Table 1). A A—A 
or an I — K lower cut-off of 5 was used to select these 
EROs. Given the heterogenous nature of the fields from 
which these objects were taken, the sample as a whole 
is not complete or uniform. However, a subset of the ob¬ 
jects observed, the EES sub-sample (Elston et al. 2001) 
does form a complete sample of 14 EROs, selected such 
that A < 20 and / — A > 5, and are listed separately 
in Table 1. Hereafter, the term ERO will imply all 27 ob¬ 
jects listed in the table and the EES sample will include 
only the 14 objects forming a complete sample. The av¬ 
erage R — K and I — K colours of the EES sample and 
the rest of the EROs are statistically indistinguishable. A 
sub-sample of the EROs selected were observed at 1250 
/mi using the IRAM 30 m telescope (15 objects) and at 
850 /(in (21 objects) and 450 /tm (5 objects) using the 
SCUBA at the JCMT. Some of these objects have also 
been observed either at 4.8 GHz or at 1.4 GHz (11 and 7 
objects respectively) using the VLA. Due to various tele¬ 
scope scheduling constraints, the same objects could not 
be observed at all wavebands. 

2.1. NIR observations of the EES complete sub-sample 

Elston et al. (2001) have completed a BRIzJK field sur¬ 
vey (hereafter, EES survey) covering ~100 arcmin 2 over 
four areas of the sky at high galactic latitudes, down 
to A^21.5. These observations were carried out at the 


in the NIR (I — A'=6, 
Cimatti et al. (1998a) 
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4 m telescope of the Kitt Peak National Observatory. 
Optical imaging was obtained with the PFCCD/T2KB, 
which gives 0.48" pixels over a 16' field. The IR imaging 
was obtained with IRIM, in which a NICMOS3 HgCdTe 
array provides 0.6" pixels over a ~2.5' field. Details of 
the observations and data analysis are given in Elston et 
al. (2001). The BRIzJK survey images were co-aligned 
and convolved to the same effective PSF of FWHM~1.5". 
Calibrations of the optical and IR images onto the Landolt 
and CIT systems were obtained using observations of 
Landolt and UKIRT standard stars, respectively. 

A catalogue of objects within each of the four fields 
of the EES survey was obtained from the K images using 
SExtractor (Bertin & Arnout 1996). Object detection was 
performed down to a level corresponding to 1.5 a above 
the sky level, with a minimum object size of 1.2 square arc- 
seconds. Photometry was obtained through 3.3" diameter 
apertures using this catalog on the BRIzJK images; the 
4 a limit in the K band is 21.4. Details of these catalogues 
are available in Elston et al. (2001). 

2.2. SCUBA observations at the JCMT 

A sub-sample was observed with Submillimeter Common- 
User Bolometer Array (SCUBA) instrument on the James 
Clerk Maxwell Telescope (JCMT) at 450 /jm and 850 pm 
in the standard point-source photometry mode during dif¬ 
ferent observing runs from 1998 to 2000. The typical opac¬ 
ity r was in the range 0.2-0.5 and 2-4 at 850 pm and 450 
(im respectively. The on-source integration times varied 
between 900 and 3600 seconds. The data reduction was 
performed using the Starlink SURF software (Jenness & 
Lightfoot 1998). For each integration, the measurements 
in the reference beam were subtracted from those in the 
signal beam, rejecting obvious spikes. Flat-field correc¬ 
tions were applied to each observation which were sub¬ 
sequently corrected for atmospheric opacity. The residual 
sky background emission was then removed using the me¬ 
dian of the outputs of the different rings of the bolometer 
as a background estimate. The flux calibration of the data 
was performed mainly using a primary calibrator (Uranus 
or Mars), thus yielding a 10% accuracy for the flux density 
scale. A poorer calibration accuracy (20%) was obtained 
for EES Lynx 2 and SA57 2-3 for which the secondary cal¬ 
ibrators CRL618 and HL Tau were used respectively. The 
individual reduced and calibrated observations were con¬ 
catenated for each source thus obtaining a final co-added 
data set. 

2.3. I RAM 30 m telescope observations 

The 1.25 mm data reported here were taken with 
the MPIfR 37-cha nnel and 19-channel bolometers 
Kreysa et al. (1998 )| at the focus of the IRAM 30m an¬ 
tenna (Pico Veleta, Spain) during observing runs in March 
1998 and December 1998 respectively. The filter set used 
combined with atmospheric transmission produces an ef¬ 


fective wavelength of about 1.25 mm; the beam size is 
11" (FWHM) and the chop throw was set at 50" and 30", 
during the first and second run, respectively. The expected 
average sensitivity for each channel, limited principally by 
atmospheric noise, was 60 mJy/vRIV, where BW is the 
bandwidth used, in Hz. However, the 37-channel bolome¬ 
ter observations were noisier than expected and therefore 
the 19-cliannel bolometer was used for the second observ¬ 
ing run. The effect of sky noise was reduced substantially 
by exploiting the correlation between signals from the dif¬ 
ferent channels using the standard three beam (beam¬ 
switching + nodding) technique, resulting in an average 
rms la value, after 6000 seconds of integration, of 0.57 
mJy and 0.4 mJy for the 37-channel and 19-channel data 
respectively. The typical on-source integration time was 
6000 seconds, distributed over two to three nights. The 
atmospheric transmission was monitored by making fre¬ 
quent sky-dips. The average zenith opacity was 0.13 dur¬ 
ing the observations in March, with a maximum of 0.2, 
and the average opacity was ~0.3 during the December 
observing run. Absolute flux calibration was performed 
using Uranus as the primary calibrator and using Mars 
and quasars from the IRAM pointing list as secondary 
calibrators. The different calibration measurements were 
consistent at a level of 5 % for both planets. Including the 
uncertainty in the planet temperatures, the average flux 
calibration uncertainty was estimated to be 10 %. Pointing 
was checked every hour and the average accuracy achieved 
was better than 3". 

The data were reduced assuming that the target 
sources are unresolved, i.e. the source sizes at mm wave¬ 
lengths are smaller than the size of the central channel. 
The remaining 36 and 18 channels (excluding one which 
suffered a large electronic loss) were then used to derive 
a low-noise sky estimate. The average value of the sky 
brightness, computed using these outer 35 and 17 chan¬ 
nels respectively, was subtracted from the signal in the 
central channel to derive the final flux density estimates. 


2.4. VLA observations 


The Very Large Array (VLA) was used to observe a total 
of 16 EROs distributed in six different fields. 4.8 GHz radio 
continuum emission was searched for in eleven of these ob¬ 
jects in the D configuration at 4.8 GHz, chosen as a com¬ 
promise between optimising point source sensitivity using 
a low resolution array and minimizing confusion through 
the usage of a high resolution array. The correlator was 
used in the continuum mode and data were acquired in 4 
IF bands, each of 50 MHz bandwidth. The flux calibra¬ 
tion and initial phase calibration was done using standard 


techniques described in Taylor et al. (1999) and the data 
reduction was done using standard algorithms in the soft¬ 
ware AIPS. Self-calibration of only the visibility phases 
was done for a few fields using background sources in the 
primary beam. The rms noise in the final images is within 
20% of the expected thermal noise. 
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Seven EROs distributed in two fields were observed us¬ 
ing the VLA in the B configuration at 1.4 GHz. The fre¬ 
quency of observation was chosen based on the expected 
steep spectral index of the continuum emission of EROs 
(typically —0.7) and poor sensitivity of the VLA at v < 1 
GHz. Optimising between confusion and the slightly ex¬ 
tended nature of EROs (>2") led to the B configuration 
being used (the number of sources within a synthesised 
beam area in this configuration with flux densities higher 
than lcr for 8 hours of integration is 0.015, using the rela¬ 
tion given in Langston et al. 1990). These observations 
were made in the line mode using 4 IFs, each with a 
bandwidth of 25 MHz and 7 available channels (the 50 
MHz bandwidth system was not used at this frequency 
as much higher closure errors are expected: Owen, private 
communication). The dataset were lranning smoothed fol¬ 
lowing standard calibration procedures and were imaged 
using multi-band synthesis (in order to reduce the effect 
of bandwidth smearing far from the phase centre). The 
tangent-plane approximation was used to correct for the 
effects of the array non-coplanarity and strong sources 
within two primary beams of the phase centre were iso¬ 
lated by defining multiple fields and simultaneously de¬ 
convolved using the CLEAN algorithm. Self-calibration 
using sources in the primary beam was done only for the 
visibility phases and the data were flagged based on exces¬ 
sive closure errors. An rms within 20% of the theoretical 
noise was achieved for the EES-Pisces field but the noise 
in the EES-Cetus field was substantially higher due to the 
presence of strong sources near the phase centre. Further 
details are given in Table 2. 

3. Results 

Out of the seven sources observed at 1250 /iin, one source, 
SA57-1, was marginally detected at the 3cr level. The mea¬ 
sured flux densities at the NIR positions of the EROs along 
with the lcr errors are listed in Table 3 for the 450 pm, 
850 /im and the 1250 pm data. As can be seen, none of 
the other observed sources were detected at these wave¬ 
lengths. 

We have detected 1.4 GHz radio continuum emission 
from the NIR positions of two EROs : EES-Cetus 1 and 
EES-Cetus 2 (hereafter, EESC1 and EESC2). None of the 
other sources observed at either 4.8 GHz or 1.4 GHz were 
detected and the corresponding measurements at the ex¬ 
pected NIR positions and the lcr errors are given in Table 
3. The 1.4 GHz image of the EES-Cetus field with the 
two detections is shown in Fig. 1. Though there appears 
to be some extended emission in the image of EESC1, 
given the signal-to-noise ratios of the two detections, the 
two sources are essentially unresolved. We estimate the 
likelihood of these two radio sources being the counter¬ 
parts of the NIR-detected EROs as follows : The posi¬ 
tional error in the radio positions is calculated for the two 
sources using the relations given in Rieu (1969), assuming 
the error in the NIR position to be 0.5", which is the lcr 
error between the radio and the optical reference frames 
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Fig. 1. Contour plot of the 1.4 GHz continuum image of 
the EES-Cetus field. The centres of the stars mark the NIR 
positions of the EROs. The two detections, EESC1 and 
EESC2 can be seen. The CLEAN beam is 6.5"x5.3"and 
is shown in the bottom left corner of the figure. The rms 
is 75 /tJy. The contour levels are at 3.3cr, 4<r, 6 cr, 8 cr, 17cr, 
32cr & 63cr. 

(Russell et al. 1990). The value of the likelihood ratio (LR) 
as defined by de Ruiter et al. (1977) is calculated to be 
1276 for EESC1 and 43 for EESC2, using the relation in 
Langston et al. (1990) to estimate the surface density of 
objects with Si .4 ghz >0.1 mJy. This implies that the 
a posteriori probabilities of the radio detections being ac¬ 
tual counterparts are 99.92 % and 97.73 % for EESC1 and 
EESC2 respectively (assuming that the a priori probabil¬ 
ity that the ERO does have a radio counterpart is 50 %). 
Hence the radio counterpart to EESC1 seems real but the 
identification of the counterpart of EESC2 is slightly less 
secure. 

4. Statistical properties of the sample 

Though almost the entire sample is undetected at all ob¬ 
served wavelengths, the flux densities measured at the IR 
positions of the sources can be used to compute the statis¬ 
tical flux densities of the sample and hence derive stricter 
upper limits. This weighted average flux density at 4.8 
GHz is 3.3 ± 4.6 /(Jy for a sample of 10 objects observed 
(the value for the 1.4 GHz data is —15 ± 12 pjy, and is 
far less significant in terms of the implied limits on the 
star formation rate). The statistical flux density at 850 
/im for the entire sample of 21 EROs is 1.0 ± 0.4 mJy and 
the corresponding number for the 1250 /tin measurements 
is 0.18 ± 0.12 mJy for 15 sources. Clearly, the 1250 /im 
and the 4.8 GHz data do not yield a statistical detection 
of the ERO population whereas the 850 /im data yields a 
marginal 2.5cr statistical detection. A similar exercise was 
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carried out for the EES sample as well and the average 
flux densities at 4.8 GHz, 1250 fim and 850 /irn are 3.5 ± 
4.8 ^Jy, 0.18 ± 0.13 mJy and 0.3 ± 0.6 mJy for nine, seven 
and fourteen objects respectively; there is no statistical de¬ 
tection for this complete sample either. We have quoted 
the values for the 4.8 GHz and the 1250 /jm data for the 
EES sample though not all sources have been observed, as 
the observed sources probably form a random sub-sample 
of the complete sample since they were selected based on 
telescope scheduling constraints. The 2.5cr result for the 
average 850 /im flux density for the entire sample can be 
traced solely to the non-EES sample of 7 galaxies, whose 
mean signal to noise ratio is la. Therefore it is probable 
that some of these EROs not in the EES sample might be 
detectable at 850 /im if observed further. In this paper, 
we will take the 3 a upper limit of 1.2 mJy for the average 
flux density at this wavelength for the entire ERO sample. 

The sample was also divided into two sub-samples us¬ 
ing an I — K colour cut-off (various cut-off values ranging 
from 5.3 to 5.9 were tried) and it was seen that there is no 
statistically significant difference in the average flux den¬ 
sities between the two sub-samples. This result holds for 
the EES sample as well. If the observed population is di¬ 
vided into ellipticals and starbursts using the I—K versus 
J — K diagram of Pozzetti & Mannucci (2000), described 
in Sect. 5.1, the two categories of EROs do not differ in 
their average flux densities either, within errors. 


5. Ellipticals or Starbursts ? 


Though attempts to distinguish between pure starbursts 
and elliptical galaxies among the ERO population have 
met with unambiguous results only for sources bright in 
IR and optical, or for sources with detectable mm/sub-mm 
emission, a few recent studies seem to show that probably 
> 70% of the EROs are old high redshift ellipticals. The 
evidence for this is based on both individual spectroscopic 
identifications of small samples of EROs (Cimatti et al. 
1999; Liu et al. 2000) and on high resolution images, us¬ 
ing morphological information (Stiavelli & Treu 2001) and 
through fitting the de Vaucouleurs’ law to radial profiles 
(Moriondo et al. 20 00). Additionally, a recent wide-field 
survey of EROs by Daddi et al. (2000a)| has shown that 
these objects are strongly clustered in the sky and this re¬ 
sult has been confirmed in two other fields by McCarthy 
et al. (2000). The strong clustering of the ERO population 
is added evidence that the majority of EROs are indeed 
ellipticals, since ellipticals are known to be more clustered 
than spirals, and also because of the narrow range of red- 
shift, 2=1-2.5, allowed for extremely red ellipticals (Daddi 
et al. 2000b). 


5.1. Diagnostic techniques 

In this section, the inferences that can be drawn from the 
observed multi-frequency flux densities of the EROs using 
group properties of ellipticals and starbursts are investi¬ 
gated. 


It is known that the I\ magnitudes of radio-loud galax¬ 
ies are correlated with their redshifts (Lilly & Longair 
1984). It is also established that galaxies of other types 
have fainter K magnitudes than radio-loud galaxies at 
that redshift (van Breugel et al. 1999, de Breuck 2000), 
which can be used to set upper limits on the redshift of 
EROs. Given the faintness of EROs in the K band, such 
an exercise constrains them to lie at z <5. 

If we consider the upper limit to the 1.4 GHz flux den¬ 
sity of the undetected ERO sample to be 0.1 mJy, (ex¬ 
trapolated from the 4.8 GHz upper limit. This excludes 
the EROs observed at 1.4 GHz which have varying up¬ 
per limits and excludes the two detections as well) then 
for z < 5, the rest frame 1.4 GHz luminosity is less than 
2x 10 25 W/Hz for Ha=0.7 (and less than 6x 10 24 W/Hz for 
Ha = 0). From the bi-modal 1.4 GHz luminosity distribu¬ 
tion of the IRAS 2 Jy sample (fig. 15 of Yun et al. 2001), 
it is clear that our sample cannot be differentiated into 
galaxies dominated by starbursts versus those by AGNs, 
which reflects the fact that our radio data is not deep 
enough to detect weak starbursts. 

Pozzetti & Mannucci (2000) showed that old ellipticals 
and dusty starbursts occupy distinctly different areas in 
the I — K versus J — K diagram and derived the theo¬ 
retical dividing line in this plane. Due to the faintness of 
most EROs in the K band, the error bars for the colours 
are too large for this diagnostic to be used profitably. 
Additionally, HR 10, the definitive example of a dusty 
starburst ERO (Andreani et al. 2000; Dey et al. 1999; 
Cimatti et al. 1998a), lies on the dividing line. Hence, 
more sensitive NIR photometry is needed in order to use 
this method of classification. 

For our sample, based on the few IR-optical colours 
available, and the upper limits from the radio, mm and 
sub-mm observations, it is not possible to determine the 
nature of each of the EROs individually. Instead, we now 
discuss the statistical properties of EROs, for the dusty 
starburst galaxy and the old elliptical galaxy components, 
seperately. 


5.2. Star formation properties of EROs 

The average 850 pm flux density limit can be used to de¬ 
rive upper bounds on the average star formation proper¬ 
ties of the sample. Assuming that all EROs in the sample 
are at a redshift of 1.5 (i.e., roughly at the redshift of HR 
10 ), we derive the average star formation rate to be less 
than 150 Mq yr _1 (from the relation given in Carilli & 
Yun 1999, and assuming the dust emissivity /3= 1.5) and 
the average FIR luminosity to be less than 1.6xl0 12 L 0 
(Ha=0). For a dust temperature of 20 K, the correspond¬ 
ing upper limit on the average dust mass is 6 x10 s M 0 
and for 50 K, the dust mass is less than 1.4xl0 8 Mq. For 
B a =0.7, SFR < 380 Mq yr' 1 , L Fm < 4xl0 12 L e and 
Mdust < 2x 10 9 Mq for Td us t=20 K and < 3.6x10 s M 0 
for Tdust=50 K (for a dust emissivity f3= 1.5 and mass ab- 
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sorption coefficient of 0.15 m 2 kg 1 at 800 /im; Hughes et 
al. 1997). 


However, we can assume that a fraction x of the EROs 
are dusty ULIGs which resemble HR 10 in their properties 
(Lfir ~4x10 12 Lq at z~1.5) and the rest are ellipticals 
with no sub-mm and mm emission. Then, from the esti¬ 
mated 5cr average flux densities of the ERO sample at 850 
fj ,m and 1250 pm and the measured flux densities of HR 
10 (Sgso M m=5.5 ± 0.6 mJy, the weighted average of the 
values quoted in Cimatti et al. 1998a and Dey et al. 1999 
and S 1250 /i m =4.9 ± 0.7 mJy, Cimatti et al. 1998a), the 
value of x can be computed. Such an exercise gives x < (36 
± 11) % and x < (12 ± 4) % for the 850 pm and 1250 pm 
data respectively. This estimate, though approximate, is 
consistent with other independent estimates of the ULIG 
fraction in EROs : <30 % (Cimatti et al. 1999; Moriondo 
et al. 2000). 

If we assume, as a conservative estimate, that the star- 
burst fraction in the ERO population is as much as 30 %, 
then from the surface density of EROs with R — K s > 5 
and K s < 19.2 estimated by Daddi et al. (2000a), the 
surface density of s tarbursts would be less than 725 
± 33 objects deg -2 . Smail et al. (1999a) calculated that 


the surface density of SCUBA sources with a 850 pm 
flux density greater than 0.5 mJy (which corresponds to 
the cut-off estimated in order to fully explain the Far 
Infrared Background or FIRB) is 17000 objects deg -2 
(for Sgso > 2 mJy which would explain half the ob¬ 
served FIRB, the surface density is 3700 objects deg -2 ). 
Hence the maximum overlap between R — K s > 5 and 
K > 19.2 EROs which are starbursts, and the high red- 
shift star forming SCUBA sources is 4 % (for a 0.5 mJy 
cut-off for Sgso /jin', it is less than 20 % for a 2 mJy cut-off). 


Given the lack of detectable sub-mm emission from the 
sample of objects, and the estimate of the corresponding 
fraction of HR 10-like ULIGs, the present study clearly 
shows that the ULIGs like HR 10 are rather unique objects 
among EROs and hence dusty strong starbursts are not 
the dominant component of this population. 


5.3. EROs as elliptical galaxies 

It can be seen from the upper limits to the radio lumi¬ 
nosities of the EROs (derived in section 5.1) that there 
are no radio-loud ellipticals in the sample. Hence these 
EROs must either be centre-brightened radio galaxies (or 
FR I; see Fanaroff & Riley 1974 and Ledlow & Owen 1996 
for definitions. The optical luminosities have been derived 
assuming a median redshift of 1.5 and m# derived from 
Table 1) or radio-quiet ellipticals. There are two detec¬ 
tions at 1.4 GHz, and the 4.8 GHz upper limits of 11 EROs 
scaled to 1.4 GHz are ~0.1 mJy. Assuming that >70% of 
these 13 EROs are elliptical galaxies, the detection rate of 
ellipticals for a 0.1 mJy cut-off at 1.4 GHz is calculated to 
be <22 ± 16 %. If the three EES-Cetus sources undetected 
in the radio are also included, the detection rate becomes 
<27+ 19 %. For a redshift range of 1-3, 0.1 mJy at 1.4 


GHz corresponds to a 1.4 GHz rest-frame radio luminos¬ 
ity of 2x 10 23 -2x 10 24 W/Hz (Ha=0). Also, the rest-frame 
R band magnitude for the sample is between —18 and 
—25.5 (assuming a RT-correction K(z)=1.122z for ellipti¬ 
cal galaxies, see Ledlow & Owen 1996). It should be noted 
that for the same set of parameters, using the data pub¬ 
lished in Gavazzi & Boselli (1999), Ho (1999), Ledlow & 
Owen (1996) and Auriemma et al. (1977), we derive the 
detection rate to be 1% to <3% for low redshift ellipticals. 


6. Continuum detections - EESC1, EESC2 and 
SA 57-1 

The two EROs, EESC1 and EESC2, with 1.4 GHz contin¬ 
uum detections (though the radio identification is tenta¬ 
tive for EESC2) have only upper limits to their sub-mm 
flux densities. The optical-IR SED of these objects cannot 
be used to classify them as either elliptical or starforming 
galaxies unambiguously. Further, the upper limit to the 
850 pm-1.4 GHz spectral index af 53 is +0.47 and +0.56 
for EESC1 and EESC2 respectively. 

The ratio of radio to sub-mm flux density of a pure 
starburst decreases with increasing redshift, i.e., for a 
given sub-mm flux density, the z=0 galaxies have the max¬ 
imum radio flux density. Hence for a given upper limit to 
the sub-mm flux density for a galaxy at an arbitrary red¬ 
shift, the z=0 radio-FIR correlation will yield the maxi¬ 
mum possible radio flux density. Therefore a sufficiently 
radio-bright AGN might possibly have a radio continuum 
strength higher than this value and will be easily classi¬ 
fied as a radio-loud AGN. Since that is not so for either of 
the two radio detections, they cannot be unambiguously 
classified as either an old elliptical or a dusty starburst. 
The properties of these galaxies are derived below assum¬ 
ing either of the two possibilities - that they are pure 
starbursts, or they are old ellipticals. 


6.1. EESC1 and EESC2 : 1.4 GHz detection 


If EESC1 and EESC2 are assumed to be dominated by 
star formation, then from the derived upper limits to the 
spectral index ctf 5 4 for EESC1 and EESC2, we can de¬ 
rive the upper limits to the redshifts of these objects. 
|Carilli fc Yun (2000)| have tabulated the values of redshift 
for a given value of the spectral index af 53 and have also 
tabulated the values for the ±lcr curves of redshift ver¬ 
sus a ® 5 4 . Using the derived value of a+Aa, where Aa 
is the lcr error on the value of the spectral index a, and 
the z - curve of parilli fe Yun (2000) , updated in Carilli 
(private communication)!)], we estimate the upper limits 
to the redshifts (see Carilli & Yun 2000 for details). 

The limiting redshifts of EESC1 and EESC2, derived 
as described above, are z < 1.5 and x < 2.0 respectively. 
Their radio continuum flux densities, if attributed solely 
to star formation, imply a star formation rate of about 


The updated values can be found at the URL: 


http://www.aoc.nrao.edu/~ccarilli/alphaz.shtml 
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1100 Mg yr -1 and 1600 Mg yr -1 respectively for the 
maximum redshifts derived above (these are for Ha=0.7. 
The values are about 700 Mg yr -1 and 1000 Mg yr -1 
respectively, for Oa=0), and less than this value for lower 
values of redshift (from the relation given in Carilli & Yun 
1999). Also, from the observed radio flux density, the FIR 
luminosity is calculated to be less than ~10 13 Lq , for 
^ < Zmax, for the two galaxies. 

If the two galaxies are assumed to be old ellipticals 
instead, then using the pure luminosity evolution models, 
the extreme R—K colours of the two galaxies ( R—K >6.1) 
imply that they lie at redshifts greater than 1.3 (Daddi et 
al. 2000b). For 2 between 1.3 and 5, the 1.4 GHz rest- 
frame luminosities of EESC1 and EESC2 correspond to 
the luminosity range of FR I galaxies or radio-quiet ellip¬ 
ticals. 

6.2. SA57-1 : 1.25 mm IRAM detection 7 

SA57-1 is detected at 1.25 mm with the IRAM 30m tele¬ 
scope at the 3cr level: Si. 25 mm =1.45±0.45 mJy. Assuming 
that this detection is real, it is not obvious whether this 
source is a pure starburst or whether it also harbours an 
AGN. From the 1.25 mm measurement and the upper limit 
to the radio emission, a lower limit of 0.9 can be derived 
for the redshift of the object (using the z + curve in Carilli 
& Yun 2000), as described in section 5.1. This is also the 
lower limit to z if SA57-1 is an elliptical galaxy. Using 
the K — z relation for radio-loud galaxies, the derived up¬ 
per limit to the redshift is 3.5. If this object were a pure 
starburst, the implied SFR is between 45-800 Mq yr -1 
for Ua=0 (and 70-1300 Mq yr -1 for fl A =0.7). The cor¬ 
responding Lfir is greater than 4.5 xlO 10 Lq . 

7 . Conclusions 

Motivated by the discovery of HR 10 and its star forma¬ 
tion properties, a sample of EROs was observed in order 
to detect radio, mm and sub-mm continuum emission and 
constrain the redshifts and star formation rates of ULIGs 
in the sample. One ERO was detected at 1.4 GHz and a 
possible radio counterpart was identified for another ERO 
at the same frequency. A third was tentatively identified 
at 1250 /an. Their redshifts and star forming properties 
were constrained using their radio-FIR spectral index but 
their nature could not be unambigously determined. Since 
the sources are faint, standard techniques to classify the 
sources in the ERO sample individually as ellipticals or 
starburst are inadequate. Weighted average flux densities 
were computed for the sample using measurements at the 
IR positions of the EROs and these values are 1.0 ± 0.4 
mJy at 850 /mi, 0.18 ± 0.12 mJy at 1250 /(m and 3.3 ± 
4.6 /(Jy at 4.8 GHz. We find no difference within errors 
in the weighted average values between EROs with NIR 
colours above and below an assumed / — K colour cut¬ 
off. If the sample is divided into ellipticals and starbursts 
based on the NIR two-colour diagnostic diagram, no dif¬ 
ferences in their average flux densities are seen between 


these groups. From the lack of detection of sub-mm emis¬ 
sion from any of the EROs in our sample, it is now clear 
that dusty strong starbursts, or high redshift ULIGs, are 
not the dominant component of this population. If it is 
assumed that such an ULIG population would resemble 
HR 10 in their SED properties, then such galaxies cannot 
constitute more than about 35 % of the EROs. From the 
observed source counts of SCUBA sources and the ERO 
surface density, we suggest that EROs contribute negligi¬ 
bly to both high redshift star formation as well as to the 
FIRB. On the other hand, if all EROs have similar star 
formation properties, the average dust mass is calculated 
to be less than 2 x 10 9 Mq (for Ua= 0.7 and Td U st=20 K) 
and the average FIR luminosity, less than 4 x 10 12 Lq. 
If not more than a third of the EROs in the sample are 
assumed to be starbursts and the rest are assumed to be 
ellipticals at a median redshift of 1.5, we calculate the de¬ 
tection rate of ellipticals for a 0.1 mJy cut-off at 1.4 GHz 
to be less than 22 ± 16 %. The corresponding number 
estimated in the local universe is <3 %. 

Therefore it is clear that the dominant population of 
EROs are not ULIGs similar to HR 10 but are proba¬ 
bly old ellipticals and weaker starbursts. The determina¬ 
tion of the properties and redshifts of the elliptical galaxy 
component of the ERO population is extremely impor¬ 
tant for constraining structure formation models (Daddi 
et al. 2000b). Towards this end, the upper limits to the 
radio and the sub-mm flux densities derived in this study 
should be used to plan future observations to detect these 
objects at these wavebands. 

Acknowledgements. The JCMT is operated by the Joint 
Astronomy Centre, on behalf of the UK Particle Physics and 
Astronomy Research Council, the Netherlands Organization 
for Scientific Research (NWO) and the Canadian research 
council. The VLA is a facility of the National Radio Astronomy 
Observatory (NRAO), which is operated by Associated 
Universities, Inc. under a cooperative agreement with the 
National Science Foundation. This work was based on observa¬ 
tions carried out with the IRAM 30 m telescope. IRAM is sup¬ 
ported by INSU/CNRS (France), MPG (Germany) and IGN 
(Spain). This research has made use of NASA’s Astrophysics 
Data System Abstract Service. We thank Wolfgang Tschager 
for help with the IRAM observations. Part of this work was 
done by NRM during his visit at Leiden Observatory and 
he thanks the Netherlands Research School for Astronomy 
(NOVA), Leids Kerkhoven-Bosscha Fonds and Leiden observa¬ 
tory for financial support and the staff of Leiden Observatory 
for their hospitality and stimulating discussions. NRM also 
acknowledges useful discussions with Frazer Owen and K.S. 
Dwarakanath about data analysis and K.R. Anantharamaiah 
for encouragement in moving on to higher redshift realms. 

References 

Andreani, P., Cimatti, A., Rottgering, H. .1. A., & Tilanus, R. 

1999, Ap&SS, 266, 267 

Andreani, P., Cimatti, A., Loinard, L., & Rottgering, H. 2000, 

A&A, 354, LI 

Auriemma, C., Perola, G.C., Ekers, R., et al. 1977, A&A, 57, 

41 



N. R. Mohan et al: EROs : radio and sub-mm continuum from EROs 


Bertin, E., & Arnouts, S. 1996, A&AS, 117, 393 
Blain, A. W. 1999, MNRAS, 309, 955 
Carilli, C.L., & Yun, M.S. 1999, ApJ, 513, L13 
Carilli, C.L., & Yun, M.S. 2000, ApJ, 530, 618 
Chapman, S. C., McCarthy, P. J., & Persson, S. E. 2000, AJ, 
120, 1612 

Cimatti, A., Andreani, P., Rottgering, H., & Tilanus, R. 1998a, 
Nat, 392, 895 

Cimatti, A., Andreani, P., Rottgering, H. J. A., & Tilanus, R. 
1998b, in proceedin gs of the Xth Renc ontres de Blois, “The 
Birth of Galaxies” (rstro-ph/9810001) 

Cimatti, A. et al. 1999, A&A, 352, L45 

Cimatti, A. 2000, in proc eedings of the ESO /ECF/STScI work¬ 
shop “Deep Fields” (astro-ph/0012057) 

Condon, J. J. 1992, ARA&A, 30, 575 

Daddi, E., Cimatti, A., Pozzetti, L., et al. 2000a, A&A, 361, 
535 


Daddi, E., Cimatti, A., Renzini, A. 2000b, A&A, 362, L45 
De Breuck C. 2000, Ph.D. Thesis 

de Ruiter, H.R., Arp, H.C., & Willis, A.G. 1977, A&AS, 28, 
211 


Dey, A, Graham, J.R., Ivison, R.J., et al. 1999, ApJ, 519, 610 
Elston, R., Rieke, G. H., & Rieke, M. J. 1988, ApJ, 331, L77 
Elston, R. et al. 2001, in preparation 
Fanaroff, B. L., & Riley, J. M. 1974, MNRAS, 167, 31 
Gavazzi, G. & Boselli, A. 1999, A&A, 343, 86 
Gear, W. K. et al. 2001, MNRAS, in press 
Graham, J. R., & Dey, A. 1996, ApJ, 471, 720 
Ho, L.C. 1999, ApJ, 510, 631 
Hu, E. M., & Ridgway, S. E. 1994, AJ, 107, 1303 
Hughes, D. H., Dunlop, J. S., & Rawlings, S. 1997, MNRAS, 
289, 766 

Jenness, T., & Lightfoot, J. F. 1998, ADASS VII, Albrecht R., 
Hook R. N., Bushouse H. A. (eds), ASP Conf. Series, Vol. 
145, 216 

Kreysa, E. et al. 1998, in Advanced Technology MMW, Radio, 
and Terahertz Telescopes, ed. Phillips, T. G., proceedings 
of SPIE Vol. 3357, 319, 

Langston, G. I., Conner, S. R., Heflin, M. B., Lehar, J., & 
Burke, B. F. 1990, ApJ, 353, 34 
Ledlow, M.J., & Owen, F.N. 1996, AJ, 112, 9 
Lilly, S. J., & Longair, M. S. 1984, MNRAS, 211, 833 
Liu, M. C., Dey, A., Graham, J. R., Bundy, K. A., Steidel, C. 

C., Adelberger, K., & Dickinson, M. E. 2000, AJ, 119, 2556 
McCarthy, P. J., Persson, S. E., & West, S. C. 1992, ApJ, 386, 
52 


McCarthy, P. J. et al. 2000, in proceedings o f the 
ESO/ECF/STScI workshop “Deep Fields” (|astro- 

ph/ |5oTl499| ) 

Moriondo, G., Cimatti, A., & Daddi, E. 2000, A&A, 364, 26 
Pozzetti, L., & Mannucci, F. 2000, MNRAS, 317, L17 
Rieu, N. -Q. 1969, A&A, 1, 128 

Russell, J. L., Lasker, B. M., McLean, B. J., Sturch, C. R., & 
Jenkner, H. 1990, AJ, 99, 2059 

Smail, I., Ivison, R., Blain, A., & Kneib, J. -P. 1999a, Ap&SS, 
266, 279 

Smail, I., Ivison, R. J., Kneib, J. -P., et al. 1999b, MNRAS, 
308, 1061 

Soifer, B. T., Matthews, K., Neugebauer, G., et al. 1999, AJ, 
118, 2065 

Spinrad, H., Dey, A., Stern, D., et al. 1997, ApJ, 484, 581 
Stiavelli, M., & Treu, T. 2001, in proceedings of “Galaxy Disks 
and Dis k Galaxies”, ASP conf. series 230, ed. Funes and 
Corsini (rstro-ph/OOlOlOC) 


Taylor, G. B., Carilli, C. L., & Perley, R. A. 1999, “Synthesis 
Imaging in Radio Astronomy II”, Sixth NRAO/NMIMT 
Synthesis Imaging Summer School, ASP Conference Series, 
Vol. 180 

Thompson, D. et al. 1999, ApJ, 523, 100 

van Breugel, W., De Breuck, C., Stanford, S.A., et al. 1999, 
ApJ, 518, L61 

Walsh, D.E.P., Knapp,G.R., Wrobel, J.M., & Kim, D.W. 1989, 
ApJ, 337, 209 

Yan, L., McCarthy, P. J., Weymann, R. J., et al. 2000, AJ, 120, 
575 

Yun, M.S., Reddy, N.A., & Condon, J.J. 2001, ApJ, 554, 803 







N. R. Mohan et al: EROs : radio and sub-mm continuum from EROs 


9 


Table 1. NIR and optical broadband magnitudes of EROs 


Galaxy 

RA 

Dec 

K a 

J-K 

I -I< 

R-K 

Ref 


(J2000.0) 

(J2000.0) 

(mag) 

(mag) 

(mag) 

(mag) 


EES Sample 

EES-Lynx 1 

08 h 48 m 30+82 

+44°52'51" 

19.8 

2.4 

5.4 

> 6.2 

1 

EES-Lynx 2 

08 ,l 48 m 43 3 .60 

+44°53'43" 

20.0 

2.9 

4.7 

> 6.2 

1 

EES-Lynx 3 

08' l 48 m 42+81 

+44° 54'34" 

19.8 

2.9 

5.4 

>6.3 

1 

EES-Lynx 4 

08 ,l 48 m 48 3 .93 

+44°57'09" 

18.9 

2.4 

5.4 

6.5 

1 

SA57 1 

13 ft 09 m 15+96 

+29°16'50" 

18.8 

2.1 

5.5 

6.4 

1 

SA57 2 

13 h 09 rr T9+21 

+29°20'20" 

19.4 

2.5 

5.3 

6.3 

1 

SA57 3 

13 ,l 09 m 19+80 

+29°21'05" 

18.9 

2.2 

5.1 

6.7 

1 

EES-Pisces 1 

23 ,l 12 m 52+93 

+00°57'06" 

19.8 

3.5 

6.0 

>7.4 

1 

EES-Pisces 2 

23 ^ 12 ^ 48+31 

+01°00'09" 

20.1 

1.7 

6.7 

> 6.6 

1 

EES-Cetus 1 

03' l 01 m 26+85 

+00°21'50" 

19.9 

2.3 

5.1 

6.5 

1 

EES-Cetus 2 

03^01^41+46 

+00°22'32" 

20.0 

2.8 

5.6 

5.5 

1 

EES-Cetus 3 

03^01^24+74 

+00°22'31" 

20.2 

2.3 

5.8 

> 6.0 

1 

EES-Cetus 4 

03 ,l 01 m 25+44 

+00°23'18" 

20.1 

2.4 

>6.4 

> 6.0 

1 

EES-Cetus 5 

03^01^24+97 

+00°25'28" 

20.0 

2.2 

5.4 

6.2 

1 

Other EROs from literature 

ERO-1 

13" 12 m 14+04 

+42°43'55" 

19.3 


6.2 


2 

ERO-2 

16"44 m 57+08 

+46°26'03" 

18.7 

2.7 

6.2 


3 

ERO-3 

17^21^46+14 

+50°02'51" 

19.6 


6.4 

8.4 

8 

ERO-5 

21^07^15+29 

+23°31'20" 

19.7 


> 6.6 


4 

ERO -6 

2l"07 m 15+46 

+23°31'38" 

19.8 


>6.5 


4 

ERO-7 

22 ,l 17 m 33+24 

+00°16'04" 

19.6 


> 6.0 


5 

ERO -8 

22 ^ 17 ^ 33+24 

+00°14'28" 

20.3 


> 6.0 


5 

J0905+3408 

09' l 05' m 30+54 

+34°08'09" 

18.3 



6.1 

6 

J1046—0016 

10 h 46 m 05+90 

-00°16'46" 

19.4 


5.1 


7 

J1721+5002 

17 ^ 21 ^ 44+50 

+50°02'06" 

20.0 


6.2 

6.3 

8 

J1019+0534 

10 h 19 m 31+87 

+05° 34'36" 

18.3 



>7.0 

9 

J1205-0743 

12"05 m 20+93 

—07°43'51" 

19.2 

2.2 

4.7 

5.6 

10 

J1249+3350 

12"49 m 43+17 

+33°50'04" 

19.7 


>4.9 

>6.3 

11 


“ Though the sample was chosen with K < 20, a revision of the photometry after sample 
selection resulted in a few objects with K > 20. 

1: Elston, R. et al. 2001, in preparation; 2: Cowie, L. L. et al. 1996, AJ, 112, 839; 3: Hu, 
E. M. & Ridgway, S. E. 1994, AJ, 107, 1303; 4: Knopp, G. P. & Chambers, K. C. 1997, 
ApJS, 109, 367; 5: Cowie, L. L. et al. 1994, ApJ, 434, 114; 6 : Eisenhardt, P. & Dickinson, 
M. 1992, ApJ, 399, L47 ; 7: Moustakas, L. A. et al. 1997, ApJ, 475, 445; 8 : McLeod, B. A. 
et al. 1995, ApJS, 96, 117; 9: Dey, A. et al. 1995, ApJ, 440, 515; 10: Giallongo, E. private 
communication; 11: Soifer, B. T. et al. 1994, ApJ, 420, LI 


Table 2. Log of VLA Observations 


Field 

#of 

Phase 

Date of 

rms 

CLEAN 


EROs 

Calibrator 

Obs. 

(^Jy) 

beam 

4.8 GHz VLA D Array Observations 

ERO-1 

1 

1327+434 

28 May 1999 

18 

34"xl7" 

ERO-5 

1 

2115+295 

28 May 1999 

15 

20"xl7" 

EES-Lynx 

4 

0818+423 

28 May 1999 

14 

19"xl6" 

SA57 

3 

1310+323 

29 May 1999 

14 

27"xl5" 

EES-Pisces 

2 

2323-032 

28 May 1999 

16 

25" x 16" 

1.4 GHz VLA B Array Observations 

EES-Pisces 

2 

2323-032 

12 Nov 1999 

18 

7"x5" 

EES-Cetus 

5 

0323+055 

6 Nov 1999 

75 

6.5" x 5.3" 
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Table 3. Radio, mm and sub-mm flux densities of EROs 


Galaxy 

1250/mi 

(mJy) 

error 

(mJy) 

850/im 

(mJy) 

error 

(mJy) 

450/.im 

(mJy) 

error 

(mJy) 

4.8 GHz 

(+Jy) 

error 

(+Jy) 

1.4 GHz 
(f‘Jy) 

error 

(dJy) 

EES Sample 

EES-Lynx 1 



-1.6 

1.8 



19.0 

14.0 



EES-Lynx 2 



-0.1 

2.0 



10.6 

14.0 



EES-Lynx 3 



1.2 

1.7 



5.4 

14.0 



EES-Lynx 4 

-0.1 

0.3 

-0.8 

2.6 



-4.9 

14.0 



SA57 1 

1.5° 

0.5 

2.9 

1.6 



-12.0 

14.0 



SA57 2 

0.8 

0.4 

3.5 

2.3 



-5.0 

14.0 



SA57 3 

-0.3 

0.4 

-5.6 

2.7 



16.0 

14.0 



EES-Pisces 1 



-1.9 

1.9 



2.6 

16.0 

-11.0 

18.0 

EES-Pisces 2 

-0.3 

0.2 

-1.0 

1.9 



-1.5 

16.0 

-27.0 

18.0 

EES-Cetus 1 



0.7 

2.4 

11 

45 



520.0“ 

75.0 

EES-Cetus 2 

0.2 

0.4 

2.4 

2.8 

128 

118 



380.0“ 

75.0 

EES-Cetus 3 



-0.9 

2.1 

5 

29 



-14.0 

75.0 

EES-Cetus 4 

0.5 

0.4 

1.9 

2.1 

17 

23 



70.0 

75.0 

EES-Cetus 5 



1.8 

2.3 

33 

32 



47.0 

75.0 

Other EROs from literature 

ERO-1 

0.7 

1.4 

1.7 

1.1 



0.4 

18.0 

<14“ 


ERO-2 

-0.8 

1.3 

1.1 

4.7 







ERO-3 

0.1 

1.1 

-0.9 

3.9 







ERO-5 



2.2 

1.1 



<70 c 




ERO-6 



3.1 

2.9 







ER.O-7 



4.2 

2.2 







ER.O-8 



-0.5 

2.1 







J0905+3408 

0.1 

0.9 









J1046—0016 

0.4 

1.7 









J1721+5002 

0.4 

1.1 









J1019+0534 

-1.6 

1.7 









J1205—0743 

-0.7 

1.7 









J1249+3350 

-0.7 

1.3 










° The three entries in bold face are detections. The rest are the flux densities measured at the NIR positions 
of the EROs, with the associated la error. These values have been used to calculate the statistical flux 
densities of the sample (see Sect. 4). 
b 3<j upper limit : Richards, private communication. 

c The upper limit to the 4.8 GHz flux density of ERO-5 is quoted as this source is 1.5 synthesized beams 
away from a strong confusing source. 





































